Exact solution of the hydrodynamical Riemann problem with nonzero tangential velocities and the 
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We give a solution of the Riemann problem in relativistic liydrodynamics in the case of ultrarelativistic equa- 
tion of state and nonvanishing components of the velocity tangent to the initial discontinuity. Simplicity of the 
ultra-relativistic equation of state (the pressure being directly proportional to the energy density) allows us to 
express this solution in analytical terms. The result can be used both to construct and test numerical schemes 
for relativistic Euler equations in (3 + 1) dimensions. 
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I. INTRODUCTION 

Solutions of the Riemann problem in the relativistic hy- 
drodynamics are of crucial importance for the construction of 
modern numerical schemes designed to solve relativistic Eu- 
ler equations. In most such schemes, it is the Riemann solver 
(usually an approximate one) that is responsible for the accu- 
racy of the method and a proper resolution of possible shock 
waves yj. Moreover, quite recently Aloy and Rezzolla ap- 
plied the analysis of solutions of the relativistic hydrodynami- 
cal Riemann problem to explain a boosting mechanism occur- 
ring in astrophysical jets, proving that the importance of such 
solutions is not merely academic |2]. 

Here by the Riemann problem we understand a Cauchy 
problem for a hyperbolic system of partial differential equa- 
tions, where initial data consist of two constant states sepa- 
rated by a discontinuity in a form of a plane surface. In the 
case of hydrodynamics, such an initial discontinuity decays, 
giving rise to three possible elementary waves: a shock wave, 
a rarefaction wave, and the so-called contact discontinuity. 
The solution of the Riemann problem is thus a non-trivial one, 
and its precise form requires investigation. 

The relativistic shock-tube problem, i.e., a Riemann prob- 
lem with zero initial velocities, was investigated by Thomp- 
son in 1 3]. Later, the Riemann problem in one spatial dimen- 
sion was solved for the ultra-relativistic equation of state by 
Smoller and Temple |4] and for the perfect gas equation of 
state by Marti and Miiller |5]. The latter work was general- 
ized by Pons, Marti and Miiller to the case, in which the fluid 
is allowed to move in the direction tangent to the discontinuity 
10], but, due to the complexity of equations, the solution had 
to be computed numerically. Then, in |7-9] Rezzolla, Zanotti 
and Pons introduced a particularly convenient way of classify- 
ing the solutions, based on the relative velocity between both 
Riemann states. 

An important progress was made by Giacomazzo and Rez- 
zolla, who analyzed the Riemann problem in relativistic mag- 
netohydrodynamics 1,10,1 . A numerical code written for this 
analysis was then used to obtain some test solutions of the 
Riemann problem for equations of state other than that of per- 
fect gas [11]. 

In this paper we present an analytic solution for the Rie- 
mann problem with non-zero velocities tangent to the initial 
discontinuity and the ultra-relativistic equation of state. This 
common equation of state is exceptional in the sense that it 
cannot be expressed in terms of the baryonic (rest mass) den- 



sity and the specific internal energy, but it relates the pressure 
directly to the energy density. This fact prevented us from a 
straightforward application of the existing numerical schemes 
solving relativistic Riemann problem. As a benefit we got 
a solution that can be expressed almost entirely in analytical 
terms. 

Solutions of the Riemann problem in which tangential ve- 
locities do not vanish can be used to construct general nu- 
merical schemes that solve equations of hydrodynamics in all 
three spatial dimensions. A solver of this kind has been im- 
plemented using the solution discussed in although it is 
not given in analytical terms, and, in order to obtain such a 
solution, one has to integrate a certain ordinary differential 
equation numerically. In the case presented here, the appro- 
priate ordinary differential equation was solved analytically, 
so the implementation of the exact Riemann solver is straight- 
forward. 

We should also note that the effects caused by the pres- 
ence of the tangential velocities in the Riemann problem are 
purely relativistic. In Newtonian hydrodynamics they do not 
influence the behavior of the solution in the direction normal 
to the discontinuity. Thus, in order to extend a given one- 
dimensional solution to the case with non-zero tangential ve- 
locities, it is only required to compute the values of those 
velocities in the intermediate states. In relativistic hydrody- 
namics all velocities couple together through Lorentz factors, 
and the presence of tangential velocities changes the solution 
quantitatively. 

Throughout this work we will assume that the reader has 
a basic knowledge of the Riemann problem for general sets 
of nonlinear hyperbolic equations (a good introduction can be 
found in fl2^). In Sees. II and III we will review basic equa- 
tions constituting our problem. Afterward, in Sees. IV-VI we 
will discuss the structures of rarefaction waves, shock waves 
and contact discontinuities respectively. Next, in Sec. VII the 
solutions of the Riemann problem will be presented, and in 
Sec. VIII we will compare them to the solutions obtained for 
the perfect gas equation of state. A summary of the paper will 
be given in Sec. IX. 



II. RELATIVISTIC EULER EQUATIONS AND THE 
EQUATION OF STATE 

The equations expressing the conservation of the energy 
and momentum in relativistic hydrodynamics are usually writ- 
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ten in the following compact form 

d^Ti"^ = 0, (1) 

where the energy-momentum tensor is that of perfect fluid, 
namely 

T"" ^{p + p)u''u'' + prf\ (2) 

Here p denotes the energy density, p is the pressure, are 
the components of the four-velocity of the fluid, and rf^ - 
diag(-l, -hi, +1, +1) is the metric tensor of the Minkowski 
space-time. Throughout this paper Greek indices will refer to 
space-time dimensions {p - 0, 1,2,3), while Latin ones will 
be reserved for spatial dimensions (; = 1,2,3). We will also 
work in Cartesian coordinates, where - (t, x, y, z). 

In order to solve the Riemann problem for Eqs. ([T]i it is 
convenient to rewrite them in the form where the derivatives 
with respect to time and spatial coordinates are separated ex- 
plicitly. To this end, we introduce the Lorentz factor W = m° 
and components of the three-velocity v' - u'/W. Due to the 
normalization of the four- velocity rf^vU^t/ = -1, the Lorentz 

factor can be written as W = 1 / Vl'^~'vy . In the above terms 
Eqs. ([TJ can be expressed as 

d,\i + 5,F'' = 0, (3) 

where 

U = ((p -H p)W^ -p,(p + pW^v\ ip + p)Wh^, 

(p + p)W^v'f , (4) 

and 

F'' = ((p + p)Wy, (p + p)W^v'v^ + d'^p, 

(p + p)W\'v^ + 6'^p, (p + p)wWv^ + 6'^pf . (5) 

Here 6'-' denotes the Kronecker's delta. 

By the ultra-relativistic equation of state we understand a 
relation p = c^p, where Cs e (0, 1) is a constant playing the 
role of the local speed of sound (for a photon gas or a gas of 
neutrinos = 1/3). This form of equation of state is com- 
monly used in cosmology; this is also the equation of state 
assumed in f^. For the ultra-relativistic equation of state (a 
barotropic equation of state of the form p = pip), in general) 
Eqs. ([1]) constitute a complete set of equations of hydrody- 
namics. 

On the other hand, the perfect gas equation of state, ex- 
ploited in most of numerical simulations in relativistic hydro- 
dynamics and in jslQ], has the form p - (y - l)«e, where y 
is a constant, n is the so-called baryonic (or rest mass) den- 
sity and e denotes the specific internal energy. The baryonic 
density is assumed to be a function satisfying the following 
continuity equation 

d^iinu'') = 0, (6) 

and the specific internal energy is defined as e - (p - n)/n. 
Thus, for the perfect gas equation of state the equations of 
hydrodynamics consist of Eqs. ([T]i and (|6]l. 



For some physical situations the baryonic density is much 
smaller than the energy density and p - n + ne ~ ne. In this 
case relations p - c]p and p - (y- \ )ne should be equivalent, 
provided that y - I - c]- The equations of hydrodynamics 
suitable for these two equations of state are, however, differ- 
ent and the solutions can differ even qualitatively (there is, 
for instance, no contact discontinuity for the ultra-relativistic 
equation of state and no tangential velocities in the Riemann 
problem, as the pressure is directly proportional to the energy 
density, and such a discontinuity is present in an analogous 
solution for the perfect gas equation of state, where the same 
pressure can correspond to different values of baryonic den- 
sity). A careful inspection of solutions of the Riemann prob- 
lem in both cases shows that they tend to each other in a suit- 
able sense. It should, however, be noted that in our case of 
ultrarelativistic equation of state the solution of the appropri- 
ate Riemann problem can be found analytically, whereas it 
was not possible for the case of perfect gas equation of state 
16]. 

In this paper, we specialize to the ultrarelativistic equation 
of state, although many results are more general, valid for 
barotropic equations of state p - pip). 

III. RIEMANN PROBLEM 

Without loss of generality, we will assume that the initial 
discontinuity is perpendicular to the x axis. Thus, neglecting 
the derivatives with respect to y and z, we can write equations 
for the Riemann problem as 

d, {{p + p)W^ -p) + d, {(p + p)W\'') = 0, (7) 

dr {(p + p)W^v'^) + 5, ((p + p)W\v'f + p) = 0, (8) 

d, {(fi + p)w\'y) + d, ((p + p)w\'\'>) = 0, (9) 

d, ((p + p)wV~) + [{p + p)W^v''v) = 0. (10) 

The structure of solutions of the relativistic Riemann prob- 
lem is exactly the same as in the corresponding Newtonian 
case, and it is, in fact, shared by general sets of hyperbolic 
conservation laws (cf. |fT2ll ). Let the initial discontinuity be 
located at x - 0, and let L and R refer to the left and right 
Riemann states, that is data for x < and jc > 0, respectively. 
The form of Eqs. (TTt-lfTOl) and the symmetry of the initial data 
suggest a self-similar solution depending on x and t through 
^ - xjt only. The initial state LR decays into three possi- 
ble elementary self-similar waves separated by some constant 
states. A smooth elementary wave, the so called rarefaction 
wave will be further denoted by '??^(^), where the subscript 
arrows refer to the direction from which particles of the fluid 
enter the wave. The other two elementary waves are discon- 
tinuities: a shock wave, denoted by S^(^), and a contact dis- 
continuity C. We will also use the symbol 'W^(^) to denote 
a shock wave S^{^) or a rarefaction wave '^^(^), when the 
actual character of the wave is not important. 

The decay of the initial state LR can be symbolically written 

as 

LR L-W^L.CR.'W^R, (11) 
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which corresponds to four different cases with "W 



S^(^j or 



The distinction between a contact discontinuity C and a 
shock wave S is based on the behavior of the pressure and 
the normal velocity (y' in our case) across the discontinuity. 
They are assumed to be continuous at the contact disconti- 
nuity and exhibit a jump at the shock wave. Since for the 
ultrarelativistic equation of state the pressure is directly pro- 
portional to the energy density, the only quantities that can be 
discontinuous across a contact discontinuity are the tangen- 
tial components of the velocity (v' and v'). This also means 
that the pressure, the normal velocity v ', and, in case of ultra- 
relativistic equation of state, the energy density are the same 
in both intermediate states L, and R^. 

The strategy of finding of the solution of the Riemann prob- 
lem can be now summarized as follows. We start by consider- 
ing a left moving wave and obtain the relation between 
the energy density p£, and the velocity in the region be- 
hind such a wave. Next, we repeat the same calculations for 
the right moving wave 'W^, to obtain an analogous relation 
between the energy density and the velocity . Since the 
energy density and the normal velocity are the same in both 
intermediate states, they can be computed from the equation 
pL,(yL ) ~ PrX'^r )■ The solution to this equation also iden- 
tifies the actual character of both waves and "W^ (the 
so-called wave pattern). We will discuss the actual forms of 
the relation between the energy density and the normal veloc- 
ity for all kinds of simple waves in the forthcoming sections. 



IV. RAREFACTION WAVE 

Let us first consider a rarefaction wave, that is, a smooth 
self-similar solution depending on t and x through ^ - xjt 
only. In this case Eqs. (l7])-(fT0l) reduce to 





= ^((p + 79)WV), 


(12) 


{(p + p)W\')-- 


-.±[(p + p)wHv'^f + p), 


(13) 


{{p + p)w\y)-- 




(14) 


((p + p)wV) = 


-.^[(p + p)W\'W). 


(15) 



Nontrivial solutions of these equations exist only if the 
Wronskian of the above set of equations vanishes, i.e., when ^ 
are the eigenvalues of the Jacobian dF ydV. Such eigenvalues 
can be easily found by exploiting the following observation. 
Let E = (p,v\v>\v-), y{ = and S = dWldT.. For a 

barotropic equation of state where p - p(p) and c\ - dp I dp, 
the determinant of y\. reads 



det^ = W**(p + pf{\- v,v'c^) . 



(16) 



Since it is positive for c, e (0, 1), the matrix J{ is invertible 
and = 'B:J\r\ Consider 

det ' - ^ j) det:?l = det (S - ^^) , (17) 



where I denotes the identity matrix. Since det^ + 0, it is 
clear that the values of ^ satisfying det (S - f >?() — are the 
eigenvalues of the Jacobian (9F V<9U. They can be easily com- 
puted to yield 



^(1 - c2) + C, V(l - V,V0 (1 - V,V'C? - (V^)2(l - C?)) 



1 - VjV'C^ 



(18) 

The eigenvalue is twofold degenerate. The expression for 
^± can be also written as 



±A 
1 + v^A ' 



(19) 



with A ^ = 1 + (l - (v*)^) (1 - Cj)/Cj, where we recognize 
the relativistic composition law for velocities. The values ^+ 
and ^- correspond, respectively, to the signals propagating to 
the right (towards larger values of x) and to the left with re- 
spect to the local flow of gas. For v'' = v" = we have A-Cs 
so that Cj can be identified with the local speed of sound. It is 
also worth pointing out that the same expressions can be ob- 
tained in a case where the pressure depends on the baryonic 
density n and the specific internal energy e with, c\ - dp I dp 
being replaced by 



1 lldp 



h \\ dn 



(20) 



where h - 1 + e-i-/;/nis the specific enthalpy. In the latter 
case, however, we are dealing with 5 instead of 4 equations, 
and the value - v " is threefold degenerate |6]. 

It can be deduced from Eqs. (l7])-(fT0b that, as long as we are 
interested in a smooth solution, the entropy density s defined 
by 



s - iiexp 



r dp' 
Jp, p'+p(p') 



jp, p'+p(p') 

where si and pi are constants, satisfies the equation 



and thus 



d, (sW) + d, {sWy"") = 0, 



(^W') = ^ {sWv") . 
d^ d^ 



(21) 



(22) 



(23) 



It should be pointed out that the entropy density given by 
Eq. (|2TJ is not conserved for discontinuous solutions, that is 
the Rankine-Hugoniot conditions following from Eq. (l22l i are 
not satisfied. In the case of ultrarelativistic equation of state 
the integral appearing in Eq. ( ITTT i can be evaluated to yield 
s - Cp'^'^'^'^'*, where C is a constant. Combining Eq. (fT4l) and 
(l23T l gives 



(^-v")— (pW) 



(24) 
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where k - c]l{\ + Cj). A similar result holds for Eq. ( fTsT i and 
V', so that for ^ + we obtain 

p''Wv^ - const, p^Wv' - const. (25) 

Let us introduce the tangential velocity v' as v' - 



^/(yyyTJv^. It follows that v' - aW^^p^", where a denotes 
a constant. Thus, from the definition of the Lorentz factor we 
have 

W^[l- (v'f) = 1 + a^p-^' = R(p). (26) 

A little longer calculation shows that 

(^ - v') W^dv' = (1 - ^v-') d Inp". (27) 

Inserting the expression for ^ = ^± into this equation and per- 
forming some algebra, one can arrive at the following relation 



+ T = z dlnp". (28) 



1 - (v-'Y 



Rc, 



Both sides of this equation can be integrated, but the precise 
form of the result depends on the value of the constant a. For 
a = (no tangential velocities) we obtain 



1 



1 - 



Cip 



(29) 



For non-zero tangential velocities one gets 



'1+v' 



±1 



= C2 



'l + VI +(1 -c?)aV 
1 - VI +(1 -c^Ja^p-^" 
Cs- Vl +(1 -cl)a^p- 
Cs + ^Jl +(1 - cl)a^p- 



(l/c,) 



2k 



2k 



(30) 



Knowing the state ahead the rarefaction wave we can thus 
compute the appropriate integration constant (Ci or C2) and 
obtain the solution in the region behind the front of the wave. 
The characteristics corresponding to this solution, treated as 
curves in the (f, x) space, form a "rarefaction fan," in which 
each characteristic correspond to a different value of ^+ (for 
the right moving wave) or ^_ (for the left moving one). For 
^ = ^0 - v'-*" we obtain dp/d^ = dv^/d^ = 0, so that the "fan" 
of characteristics originating at the discontinuity has a "zero 
opening angle." Remarkably, Eqs. (fT2b-(fT5b give no condi- 
tions for v'' and v'' in this case. This corresponds to the contact 
discontinuity which will be treated later in this paper. 



V. SHOCK WAVE 
Rankine-Hugoniot conditions for Eqs. ^ can be written as 

[[rn]«, = o, (31) 

where n'^ is the unit vector normal to the surface of discon- 
tinuity and [[/]] represents the jump of a given quantity / at 



the discontinuity. Since we are interested in establishing the 
state behind the wave basing on the state in front of it (the 
left or the right state in the Riemann problem depending on 
the direction in which the wave propagates), we will adopt a 
notation in which values referring to the state in front of the 
shock wave are denoted with a bar while unaltered symbols 
are reserved for the values behind the discontinuity. In such a 
notation, a jump of a quantity / reads [[/]] - f - f (a similar, 
simplified notation was used in jlSlil ). 

Assuming that the discontinuity surface is a plane normal to 
the X axis, we can write components n'^ as n'^ = Wj(ys, 1,0, 0), 
where Ws - II ^f\~-V^. The quantity Vj has a natural inter- 
pretation of the coordinate velocity of the discontinuity. In 
this case Rankine-Hugoniot conditions have the following al- 
gebraic form 

ffpH'^.^^jjy^^ [[^^2^.jj^ (32) 

[[pwV]]y, = [[pivV)' + ^pll, (33) 

[[pW^v' j] y, = [[plV^V]] , (34) 
[jpW^/jjy.^jjplV^vV]], (35) 

where we have assumed an ultrarelativistic equation of state. 

In the case of zero tangential velocity only Eqs. (|32]|-(|33]| 
are relevant. The shock wave velocity can be expressed as 

y. = [[pWV]]/[[pW2-^p]], (36) 

which, inserted into Eq. ( [33T l. gives 

[[pW^'v']]" = [[pW2(v-)2 + ^p]] [[pW2 - Kp^ . (37) 

The above equation yields 



(p/p)2-2(©+l)(p/p) + l = 0, 



(38) 



where = W^W^iy " - v^)^/(2a-(1 - a-)), and the only physical 
solution for p is given by 



p = p(i + 0+ V(i + 0)^ - 1)- 



(39) 



This equation, similarly to the rarefaction wave described 
above, gives the relation between the post-shock density p and 
the post-shock velocity v' . 

For a case with non-vanishing tangential velocity a similar 
calculation can be done. We start by multiplying both sides of 
equation (l33T l by Vs and add the result to equation (|32] |. Then, 
the expression for pW^ can be written as 



(v^ - V,)(l -y-y,)' 

where we have assumed that Vs + v'. Eqs. ( l34l l. 
following expression for the square of v' 



Ti\2 



(y, - v^Yp^w\v') 

p2W4(y, - V^)2 



(40) 
give the 

(41) 
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Inserting these two results to Eq. (bZb yields an equation 
which, after suitable rearrangement of terms, can be written 

as 
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(1 -Fy,)(i-y-y,)(v^-y,) 
x{(i-v-*y.) [(i-v'"y.)(i-v^y.) 



(v')'(i-v'^y.)(i-yf) 



(42) 



0. 



Physical values of can be now expressed in terms of v ' as 
the solutions of the cubic equation 



(1 - v^y,)(i - vv,) - i(v^ - y,)(v-'' - y,) 



(43) 



-(v^d - y-y,)(i - yf) = 0. 



This can be done, for instance, by using one of the Cardano's 
formulae. Finally, by combining Eqs. (HTt and ( l40l l we can 
obtain the following expression for the post-shock density p 
as the function of the post-shock velocity v* 

-pw\v' - y,) [(1 - {v'f){\ - v>,y,)2 - (v')2(i - v-'y,)2] 



(v»-yo(i-v"^y,)(i-v-'y.,) 



VI. CONTACT DISCONTINUITY 



(44) 



For Ys - I'' Eqs. (I32b-(l35]l have a non-trivial solution 
where - v*(= and p - p, while velocities and can 
exhibit an arbitrary jump. This corresponds to the so called 
contact discontinuity — the one co-moving with the fluid. 

Obviously, such a discontinuity can only be present in case 
of non-vanishing tangential velocities. In this respect, there 
is a qualitative difference between solutions of the Riemann 
problem for the ultrarelativistic equations of state and those 
obtained for the perfect gas equation of state. In the latter 
case only the pressure and the normal velocity v" have to be 
continuous across the contact discontinuity, and there is no 
such requirement for the baryonic density and the specific in- 
ternal energy. Thus, in case of the perfect gas equation of 
state, one usually observes the contact discontinuity also in 
a strictly one-dimensional problem (with vanishing tangential 
velocities), and such a discontinuity is absent in analogous so- 
lutions with ultra-relativistic equation of state. 



VII. SOLUTIONS OF THE RIEMANN PROBLEM 

The distinction between a shock and a rarefaction wave is 
based on the relation between the pressure in front and behind 
the wave 1 14]. If the pressure p in front of the wave is larger 
than the pressure p behind it, we are dealing with a rarefac- 
tion. The converse case with p > p corresponds to a shock 
wave. Let p = >S^(^)(y*) denote the post-shock energy den- 
sity p understood as a function of the post-shock velocity v* , 
as it can be computed from Eq. ( |44| ). As usual, the directions 



,-/ = 0,TV^ 
' = 0.5, TV^ 
' = 0.8, TV^ 
■- 0.865, TV-H 

v' = 0,'W^ 
' = 0.5,'W^ 
' = O.S,'W^ 

-- 0.865, nv^- 



-0.8 -0.6 -0.4 -0.2 







0.2 0.4 0.6 O.i 



FIG. 1. The dependence of the energy density p on the velocity 
behind the wave for the ultra-relativistic equation of state with = 
1/3. Different curves refer to values of the tangential velocity v' 
in front of the wave equal to 0, 0.5, 0.8, and 0.865. The velocity 
v' in front of the wave was equal 0.5, and the density p was set to 
1. Increasing curves correspond to the right moving waves, while 
decreasing ones to the left moving waves. 



of the arrows correspond to the direction from which the fluid 
enters the wave. A similar function giving the energy density 
behind the front of the rarefaction wave will be denoted by 
p - ??^(^)(v'). It follows from results of the preceding sec- 
tions that the general expression for the energy density behind 



a wave "W- 



can be written as 



p = nV^(v^) 



for a right moving wave, and 



(45) 



(46) 



for a left moving one. Here v" refer to the velocity in front of 
the wave. Such functions are illustrated on Fig.[T]for diff'erent 
values of the tangential velocity in front of the wave v'. 

Given two initial states L and R we can always compute 
both functions p - 'W^iy'') and p — ^^{v"), and find the 
intersection of their graphs. This occurs for some vl and p, 
common for both intermediate states L, and R^. Such an in- 
tersection has been depicted on Fig.|2]for some arbitrary states 
L and R. 

In order to complete solving the Riemann problem, one 
only has to find locations of the interfaces between different 
states in the solution. The location of the shock wave S is 
given by its speed Ys, which can be easily computed after the 
value of x'l has been established. The contact discontinuity 
C, dividing both states L, and travels with the velocity 
vj. The velocity of the head of the rarefaction wave is given 
by the expression for ^± (plus for "R^, minus for "R^) com- 
puted for the suitable Riemann state. The location of the tail 
of the rarefaction wave can be established by the condition 
that the velocity in the rarefaction wave should reach the 
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•R^(v') 
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FIG. 2. Intersection of the graphs of p(v') for left and right moving 
waves. Here both Riemann states correspond to v' = 1/2, v' = 1/2. 
The energy densities in both states differ: p/. = 10, Pr = 1. The 
curves are computed for the ultra-relativistic equation of state with 
= 1/3. 
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FIG. 3. Time snapshot of the solution of the Riemann problem for 
t = 1. The left initial state is given by = 1, = 1/2, = 1/3 
and the right state by pjt = 20, v^^ = 1/2, and vj^ = 1/2. 



value of p,. The velocity of the tail is given by ^± computed 
for the suitable intermediate state (adjacent to the rarefaction), 
but the straightforward application of formula (fTsT i requires a 
prior calculation of v' in this state. The values of v' in both 
intermediate states can be easily computed from Eq. (l40l l (for 
the state behind the shock wave) and from Eq. (|26] | (for the 
state adjacent to the rarefaction wave). 

An example of the solution of the Riemann problem for the 
ultrarelativistic equation of state with = 1/3 is shown on 
Figs. [3] and |4] Here the left initial state was given hy pi = 1, 
vl = 1 /2, = 1 /3 and the right state by pr = 20, = 1 /2, 
= 1/2. It is interesting to note the presence of a contact 
discontinuity in the tangential velocity on Fig.|4] 



FIG. 4. Time snapshot of the solution of the Riemann problem with 
the same initial data as on Fig. [3] Here the solid line corresponds to 
the velocity v', while the dotted one depicts v"^. 



Vm. COMPARISON OF THE RESULTS WITH 
SOLUTIONS FOR THE PERFECT GAS EQUATION OF 
STATE 

Solutions presented in preceding sections can be compared 
with solutions of the Riemann problem obtained for the per- 
fect gas equation of state in [5,6]. These are, in general, com- 
pletely different solutions, however, as pointed out in the sec- 
ond section, the perfect gas equation of state p - {y - l)ne 
tends to p = (y - l)p in the case where p - n + ne ~ ne. Such 
a condition can be imposed on initial data by assuming that 
n «; ne. It can be observed that the solutions of the Riemann 
problem with the perfect gas equation of state tend to those 
for the ultrarelativistic equation of state as n/e — > 0, where 
the solutions are understood as functions of the self-similarity 
variable ^ - x/t. 

It should, however, be noted that although characteris- 
tic speeds of propagation of rarefaction waves, shocks, and 
contact discontinuity tend to those obtained for the ultra- 
relativistic equation of state, they are different in each of the 
examined solutions. Thus, having a solution for the perfect 
gas equation of state which is very close to the one for the ul- 
trarelativistic equation of state in the self-similarity variable ^, 
it is always possible to consider a sufficient time f, after which 
both solutions, treated as functions of x, will vastly differ on 
an arbitrarily large subset of the domain. 

The comparison has been performed for solutions with dif- 
ferent values of initial pressures and velocities (both normal 
and tangential to the initial discontinuity) basing on numerical 
schemes provided by Marti and Miiller 1 1] and our solutions. 
In all examined cases the solutions for perfect gas equation 
of state tend to those for ultrarelativistic one in a similar way. 
An example is shown on Fig.H] where we have plotted the en- 
ergy density of a solution corresponding to the ultrarelativis- 
tic equation of state with = 1 /3 (solid line) together with 
densities computed for the perfect gas equation of state with 
y = 4/3 (dotted lines). All solutions of this example were ob- 
tained for the following initial conditions = 1 /3, = 1 /2, 
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FIG. 5. Comparison of solutions obtained for ultrarelativistic (solid 
line) and perfect gas equations of state (dotted lines). Different so- 
lutions for the perfect gas equation of state were obtained for the 
initial data corresponding to n/e equal 1.0, 0.1, 0.01, and 0.001 in 
both Riemann states. Other parameters of the initial states were set 
to PL = 1/3, vl = 1/2, v[ = 1/3, PR = 20/3, and = = 1/2. 

= 1/3, PR = 20/3, and = vj^ = 1/2 (the solution the 
for ultrarelativistic equation of state is thus the same as the 
one on Figs. [3]and|4]i. Different solutions for the perfect gas 
equation of state were computed assuming the values of n/e 
in both initial states equal to 1.0, 0.1, 0.01, and 0.001. 



IX. SUMMARY 

We have presented an exact solution of the Riemann prob- 
lem for the ultra-relativistic equation of state, with arbitrary 
initial velocities, both normal and tangential to the initial dis- 
continuity. Such a solution can be used for testing and con- 
struction of the numerical schemes which solve relativistic 
Euler equations in (3 -(- 1) dimensions. In fact, our original 
motivation for dealing with the problem presented in this pa- 
per was to provide a test solution for numerical studies of hy- 
drodynamical perturbations in the cosmology of the early uni- 
verse, where the ultrarelativistic equation of state is a frequent 
choice. 

We also point out that the boosting mechanism described in 
10] in the context of astrophysical jets is also exhibited by the 
solution of this paper. This can be observed on Fig.|4] where 
the tangential velocity v' in the region behind the rarefaction 
wave is larger than any of the velocities of the initial states. 

We have also compared our solution with a similar one ob- 
tained for the perfect gas equation of state in [5, 6] in the limit 
of vanishing baryonic density. In all examined cases solutions 
for the perfect gas equation of state, treated as functions of the 
self-similarity variable ^, tend to those for the ultra-relativistic 
equation of state. 
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